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Abstract 

We report here on the results of numerical searches for PPT states with specified ranks for density 
matrices and their partial transpose. The study includes several bipartite quantum systems of low dimen- 
sions. For a series of ranks extremal PPT states are found. The results are listed in tables and charted in 
diagrams. Comparison of the results for systems of different dimensions reveal several regularities. We 
discuss lower and upper bounds on the ranks of extremal PPT states. 

1 Introduction 

In recent years the study of entanglement in composite quantum systems has taken several different direc- 
tions. One direction is the study of entanglement from a geometrical point of view [1,2, 3, 4, 5]. This has 
led to questions concerning the relations between different convex sets of Hermitian matrices, where the 
full set of density matrices is one of them. The main motivation for the interest in these convex sets is the 
information they give about the general question of how to identify entanglement in a composite system. 
Unless the system is in a pure quantum state, the knowledge of the corresponding density matrix does not 
readily disclose the state as being entangled or non-entangled. And the complexity of the corresponding 
problem, known as the separability problem, increases rapidly with the dimensionality of the quantum 
system [6]. 

The density operators are the positive, normalized, semidefinite Hermitian operators that act on the 
Hilbert space of the quantum system, and in the following we shall use the notation V for this set. Another 
convex set is the set of non-entangled states, usually referred to as separable states, and we use S as notation 
for this subset of V. Since the set of entangled states is the complement to the set of separable states, within 
the full set of density matrices, the question of identifying the entangled states can be reformulated as the 
question of finding the boundaries of the convex set of separable states S. 

For a bipartite system there is furthermore a convex subset of the density matrices, here referred to 
as V, which is closely related to the set of separable matrices. This is the subset of density matrices that 
remain positive semidefinite under the operation of partial transposition of the matrix with respect to one 
of the two subsystems of the composite system. For short these states are called PPT states. A necessary 
condition for separability of a density matrix is that it remains positive under partial transposition, and 
thus the set of separable states is included in the set of PPT states, S C V [7]. For bipartite systems of 
dimensions 2x2 and 2x3 the two sets are in fact identical [8], but in higher dimensions the separable states 
form a proper subset of the set of PPT states. However, numerical studies have shown that for systems of 
low dimensions, like the 3x3 system, the set V is only slightly larger than S [4, 9]. 

The necessary condition that the separable density matrices remain positive under partial transposition 
is important, since this condition is easy to check. It effectively reduces the separability problem to a 
question of identifying the PPT states that are entangled, i.e., that do not belong to S. These states are also 
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interesting for a separate reason, since they are known to carry bound entanglement, which means that the 
entanglement is not available through entanglement distillation, a process where entanglement of mixed 
states is transferred to a set of pure quantum states [10]. 

All convex sets are in principle determined by their extreme points, which for T> are the pure quantum 
states and for S are the pure product states. The set V also has the pure product states as extreme points, but 
in addition there are other extreme points that are not fully known. In a previous publication [9] we have 
found a criterion for identifying extremal PPT states, and we have described an algorithm to systematically 
search for such states. By use of the method, a list of extremal PPT states were found and presented for a 
series of low-dimensional systems. 

In the present paper we follow up the study of entangled PPT states in low-dimensional systems by 
use of numerical methods. The method described here is different from that of the publication [9] in the 
sense that it does not make use of direct searches for extremal PPT states, but rather of searches for PPT 
states with specified ranks for the density matrix and for its partial transpose. By systematically searching 
through matrices of different ranks we have obtained in this way, for a series of low-dimensional systems, a 
list of low rank PPT states, many of which are identified as extremal and others as non-extremal entangled 
PPT states. These results supplement those of [9], where only a limited set of different types of extremal 
states were identified. 

The results we have obtained for different low-dimensional systems show certain regularities. In par- 
ticular we find extremal PPT states for essentially all ranks of the density operator and its partial transpose, 
when these lie between an upper and a lower limit. We suggest general expressions for these limits and 
relate them to generic properties of the image and kernel of the density matrices. A special focus is on 
the properties of the extremal PPT states of lowest rank. In a separate publication [11] we follow up these 
results by a specific study of the lowest rank extremal PPT states of the 3x3 system (which are also identi- 
cal to the lowest rang entangled states in higher dimensions [12]) . As discussed there, these states can be 
classified by a small number of parameters, and an interesting question is whether a similar classification 
can be given for other extremal PPT states . 

2 The method 

We consider a bipartite quantum system with Hilbert space H = Ha <S> where A and B label the 
two subsystems, and Ha and H b are of dimensions Na and Nb respectively. The density operators p 
satisfy the normalization condition Tr p = 1, but often it is convenient to give up this condition and rather 
consider all density operators that differ by a normalization factor to be equivalent. The set of normalized 
density operators we refer to as T>, while JC(D) is the positive cone of non-normalized density operators. 
A similar notation is used for other subsets of Hermitian operators. Partial transposition is the operation 
on density operators that corresponds to transposition of indices of one of the subsystems, pfj — > (here 
chosen as subsystem E). For density operators of the full, composite system, we refer to this operation as 
p — > p p , and it maps the convex set V into another convex set denoted V p . The operation will depend on 
the choice of subsystem ( A or E) and on the choice of basis in the corresponding Hilbert space. However, 
the distinction between these different choices is of no importance for our discussion. We note in particular 
that the mapping of sets V — > V p is independent of the choice. The set of PPT states is defined as the 
section V = VfXD F ', which means that it consists of the positive semidefinite operators that remain positive 
semidefinite under partial transposition. 

2.1 Searching for density operators of given ranks 

The Hermitian matrices define a real vector space, and it is convenient to introduce a complete set of 
matrices that are orthonormal with respect to the trace norm 

Tr(MjMj) = 5ij (1) 

A general Hermitian matrix M is then described as a vector x with real components 

Xi = Tr(MMj) (2) 
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For a composite system of Hilbert space dimension N — N^N B the vector space of Hermitian matrices is 
of dimension N' 
The algorithi 
We expand p as 



of dimension N 2 = N\N%. 

The algorithm we use to find PPT states p with specified ranks (to, n) for p and p p is the following. 



p = p(x) = X * M * ( 3 ) 

i 

The eigenvalues of p we write as Xi = Ai(x), with Af = Af (x) as the eigenvalues of the partially 
transposed matrix p p . The eigenvalues of each matrix are listed in decreasing order, and for the density 
matrix that we are searching for, a certain number of the eigenvalues should vanish. Thus we want to have 
Xk = for k = m + 1, N and X P = for k = n + 1, ...,N. The eigenvalues that should vanish we 
treat as components of a new vector /x, so that 

A* = [A m +i, A m+ 2, Aat, A n+1 , A Jl+2 , 

a£] (4) 

and the problem is then to find the point x which solves the equation jiz(x) = 0. 

We choose a starting point x such that p = p(x) as well as p p are positive semidefinite matrices (it is 
not strictly necessary that the positivity conditions hold to begin with, since they will automatically hold 
for the solution we obtain in the end). If the equation /x(x) = is not already solved, we search for a better 
approximate solution x' = x + Ax. The linear approximation to the Taylor expansion gives an equation 
for Ax, 

/i(x) + (Ax • V)/x(x) = (5) 

In matrix form the equation can be written as 

B Ax = -/x (6) 

with Bij = dpi/dxj. It implies another equation, 

AAx = b (7) 

where A = B T B is a positive, real symmetric matrix, and where b = B T fi. We use the conjugate 
gradient method to solve the last equation for Ax. The conjugate gradient method is useful because it 
works even if A is singular. Next, we replace x by x' = x + Ax, and iterate in order to get successively 
better approximate solutions. 

The matrix B is computed in each iteration by the formula 

where ipk is the eigenvector of p with eigenvalue Afe. A similar formula is used for the derivatives of 
X P . Since these formulae are valid in first order non-degenerate perturbation theory, this raises a question 
concerning convergence of the method at a point of degeneracy. However, in practice we find that the 
method works well when the dimension of the system is not too large. 

Note that, by the way the method works, all the states we find are PPT. That is the case since in the 
iterative search for a state with a certain number of vanishing eigenvalues for p and p p , the eigenvalues are 
always ordered in such a way that the lowest eigenvalues are forced to be zero. This means that both the 
density matrix and its partial transpose will be positive semidefinite. 

We have applied the method to a series of low-dimensional systems, and the results are listed in the 
tables in Appendix B and in the diagrams in the next section. The convergence of the method slows down 
with increase of the Hilbert space dimension 7Y, and in the form we have implemented the algorithm, 
the practical limit of the dimension is N < 20. For the systems we have studied, the method has been 
used repeatedly with different starting points for each choice of ranks (to, n). In most cases the iteration 
converges, but in some cases that does not happen and these iterations are then simply aborted. For some 
values of m and n we do not find any density matrix with the given ranks, and the method will then in most 
cases converge to a density matrix with lower rank for either p or p p , or both. We have not imposed any 
restriction to avoid that, so in practice the method searches for density matrices of ranks equal to or lower 
than the specified values (m, n). 
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2.2 Determining the dimension of the face of /C(P) 

For each density matrix p found in the searches, we have evaluated and listed certain properties. Among 
these are the local ranks (ta, tb) of the density matrices, defined with respect to the subsystems A and B. 
These are the ranks of the reduced density matrices pa and ps- The most interesting cases are those where 
the density matrices have full local ranks (ta, rs) — (Na,Nb). If that is not the case the density matrix 
can be viewed as belonging to a composite system of lower dimension, which is embedded in the higher 
dimensional system. 

We have also evaluated and listed the dimension of the face T of the convex cone JC{V) to which p 
belongs. The extreme points of V correspond to one-dimensional, positive rays mK{P), and are conse- 
quently characterized by dim T = 1. Therefore the extremal states can be identified in the tables as the 
density matrices with this minimal value for the dimension of the face. 

The method we use to evaluate dim J" has earlier been described in [9]. It is based on the fact that the 
face T to which p belongs can be viewed as the section between a face (or all) of JC(D) and a face (or all) 
of JC(V P ). This means that p = p(x) satisfies two equations, 

Px = x , Qx = x (9) 

with P as the orthogonal projection on the subspace in the vector space of Hermitian matrices defined by 
the face of JC(D), and Q as the projection on the subspace defined by the face of JC(T> P ). Note that P and 
Q are real and symmetric N 2 x N 2 matrices. The method for computing these matrices for a given density 
operator p is described in [9]. 

As follows from the above equations, the section between the subspaces defined by P and Q is spanned 
by the eigenvectors of the composite, real symmetric matrix PQP (or alternatively QPQ), with eigenval- 
ues equal to 1. This implies that dim F can be determined simply by diagonalizing the composite matrix 
and counting the number of such eigenvalues. 

As a consequence of the way T is constructed, there is a geometrical constraint on the possible values 
of dim T for given ranks (to, n) of p and p p [9]. Thus, the faces of K{T>) and JC(D P ) are of dimensions 
to 2 and n 2 respectively, and therefore the face of JC(D) is specified by A 2 — m 2 linear constraints and 
similarly the face of JC(V P ) is specified by A 2 — n 2 linear constraints. If these sets of constraints are 
independent they determine the dimension of the section T as dim J 7 = to 2 + n 2 — N 2 . However, if 
the two sets of constraints are not fully independent, the dimension of the section is larger. Therefore the 
following inequality is generally valid, 

dim J" > to 2 +n 2 - N 2 (10) 

For extremal PPT states, with dim T = 1, this gives an upper bound to the ranks of p and p p , 

to 2 + n 2 < N 2 + 1 (extremity) (11) 

2.3 Counting the number of product vectors 

For each density matrix p, we have examined the image (Im p) and kernel (Ker p) for the presence of 
product vectors, and the numbers of such vectors are listed. The method we use to search for product 
vectors in a given subspace of H is described in some detail in Appendix A. To briefly outline the approach, 
let us assume that we search for product vectors in Imp, with P as the orthogonal projection on this 
subspace. The vectors should satisfy the condition 

(1-P)(0®X)=O (12) 

and this can be re-expressed as a minimization problem, which we in the Appendix refer to as a double 
eigenvalue problem. The solutions of (12) are the minima of the function 

/=(0 t ®X t )(l-^)(0®x) (13) 
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with / = 0, and such minima can be found by the iterative approach described in the Appendix. By varying 
the starting point of the iteration, different minima can be identified, and by systematically searching for 
minima of / and of 

i-f = {<i>®x)*P{<i>®x) (14) 

we have reproduced and counted the product vectors in Im p and Ker p for every density matrix found in 
the numerical searches. In the tables the number of product vectors and the number of linearly independent 
product vectors are listed for Im p and Ker p. 

The same type of iterative method for minimization over product vectors has been applied in the sepa- 
rability test that is described in the next section. 

3 Discussion of the results 

3.1 Diagrammatic representation 

The results of the searches are tabulated in Appendix B, and are also included in condensed form in the text 
below, as the diagrams in Figs. 1-4. The variables of the two axes of the diagrams are the ranks to of p and 
n of p p . The small open and filled circles in the diagrams indicate the ranks for states that are found in the 
numerical searches. As discussed above, all the states produced in this way are PPT. The filled (red) circles 
represent extremal PPT states with full local ranks, while the open circles represent non-extremal PPT 
states, either separable or entangled. Note that all diagrams are symmetric under the interchange to o n, 
since in the search for PPT states there is no intrinsic difference between p and p p . 

The states we find by use of our method we refer to as typical for the chosen ranks (to, n). For some 
ranks there will exist states with untypical characteristics, which we do not find in the searches, presumably 
due to low dimensionality of the set of such states. We note in particular, when we compare the diagrams 
of the 3x3 and 4x4 systems, that most of the states of the 3x3 system are not seen in the diagram of the 4x4 
system. However, we know that all the states of the lower-dimensional system should be present, in the 
form of states with less than full local ranks. 

For most choices of ranks (to, n) the density matrices that we find in repeated searches with different 
initial conditions, are all found to have identical characteristics in the form of the parameters listed in the 
tables. There is only a small number of exceptions where density matrices with different characteristics but 
equal ranks (m, n) have been found. These are all listed in the tables. 

3.2 Different groups of PPT states 

There is a clear similarity between the diagrams of Figs. 1-4, with the states for each diagram being sepa- 
rated into groups with different characteristics. One group consists of the low-rank states with m = n < 
max{iVA, Nb}- These low-rank states, which are represented by the series of green circles in the dia- 
grams, are all separable, with equal ranks for p and p p . As shown in the tables, all the corresponding 
density matrices p contain a number of product vectors in their image which is equal to their rank to. In 
fact, it is obvious that one can construct such states by taking randomly a small number of pure product 
states and forming a convex combination of these. Such a sequence of separable states can be defined in 
all dimensions, and they are not of further interest for our discussion. 

The remaining states, with ranks to > fci&x{Na, Nb} and n > max{NA, Nb}, are restricted to the 
upper right corners of the diagrams. For these values of the ranks, all states produced in our searches, with 
one exception, have full local ranks. (The exception is the case (to, n) = (5, 5) of the 4x4 system, where in 
addition to full rank states also states with less than full rank are found.) The systems of dimension 2x2 and 
2x3, represented by the diagrams in Fig. 4, are special, since for these systems all PPT states are separable. 
However, for the higher-dimensional systems we note that the above restriction on the ranks of p and p p 
coincides with a lower bound, found by P. Horodecki et. al. [12], on the rank of entangled PPT states with 
full local ranks. (In the following we shall refer to this as the HLVC bound.) In fact, with the exception of 
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Fi gure 1: Diagrams for the ranks of PPT states in the composite systems of dimensions 3x3 and 4x4. The coordinates 
m and n on the axes are the ranks of the density matrix p and its partial transpose p p . States found by use of the 
numerical method discussed in the text are indicated by small circles or dots (filled circles). The green circles in the 
lower left corner indicate ranks for a special set of separable states, the red dots correspond to extremal PPT states and 
the red circles in the upper right corner to non-extremal PPT states. Ranks (5, 5) in the 4x4 system is an exception, 
there we find both separable states and a special type of entangled, non-extremal PPT states. The unbroken, horizontal 
and vertical straight red lines show the lower bound for entangled PPT states with full local ranks, the similar dashed 
red lines show our conjectured lower bound for extremal PPT states with full local ranks, and the red circular arcs show 
the upper bound for extremal PPT states. The dashed green 45° straight lines show the upper bound for the application 
of the separability criterion described in the text. For all ranks indicated by red circles up to and on this line (and above 
the circular arc), the corresponding states are always found to be entangled PPT states. 



some of the (5, 5) states of the 4x4 system, we find for all ranks above and sufficiently close to this lower 
bound that the states are not only entangled but extremal PPT states. 

As displayed by the diagrams, there are a few intermediate values of m, above the HLVC bound, where 
only symmetric matrices, m = n are found. We shall discuss these states separately in Sect. 3.4, and focus 
now on the states with ranks m > Na + Nb — 1 and n > Na + Nb — 1- For almost all ranks that satisfy 
these inequalities we find PPT states. The only exceptions are some of the cases with the largest asymmetry 
between the values of m and n. In fact we cannot rule out that we miss these states as a consequence of the 
method we use. We find that with large asymmetry in m and n the numerical method seems preferably to 
pick up matrices with lower and more symmetric ranks. Therefore a much larger number of searches have 
to be performed in order to find density matrices when the ranks are highly asymmetric. 

With exception for the systems 2x2 and 2x3, all the states we find with ranks between the lower bounds 
m > Na + Nb — 1 and n > Na + Nb — 1 and the upper bound m 2 + n 2 < N 2 + 1 are extremal (and hence 
entangled) PPT states. The upper bound is the constraint (1 1) on extremity that has already been discussed, 
and in the diagrams this bound is displayed as the red circular arc. For ranks above this bound, we have 
used a separability test [12] on the states, and this test shows that for all ranks that satisfy the inequality 
rn + n < 2N — Na — Nb + 2, we find entangled PPT states. Beyond this limit the test is not applicable. 

3.3 Dimensions of faces and numbers of product vectors 

The structure of the diagrams discussed above can to some extent be related to simple regularities of the 
parameters in the tables. We focus first on the list of the values of dim J 7 . This number is constrained 
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Figure 2: Diagrams for the ranks of PPT states in the composite systems of dimensions 3x4 and 3x5. States found by 
use of the numerical method discussed in the text are indicated by the small circles or dots. Also here entangled PPT 
states are found for all ranks indicated by red circles up to and on the dashed green line. For further explanation of the 
diagrams see Fig. 1. 



by the geometric bound (10), dim T > rn 2 + n 2 — N 2 . We note that for all states where the ranks are 
sufficiently large to give a positive number for this lower bound, the constraint is satisfied with equality. 
This means that the linear constraint equations (9) that define J 7 as a section between faces of JC(T>) and 
JC(D P ) are all independent. 

For all ranks m,n > Na + Nb — 1 which give negative values for the bound, we find dim T = 1. This 
is the minimum value consistent with the fact that the density operator p is located on both faces. Other 
"accidental" relations between the two intersecting faces, therefore, seem not to be present. The fact that 
dim F takes the minimal value consistent with this condition gives an explanation for why all the states 
with ranks between a lower and an upper bound are found to be extremal PPT states. For sufficiently low 
ranks, m, n < Na + Nb — 2, we find states where dim T does not take this minimum value. These states 
we therefore see as corresponding to more special constructions. 

The numbers of product vectors we find in Im p and Ker p also show a simple regularity. For a long 
sequence of high ranks m all states have no product vector in Ker p and a complete basis (in fact an 
overcomplete set) of product vectors in Im p. When lowering the rank there are in some of the lists a small 
set of intermediate ranks where there is no product vector in neither Im p nor Ker p, and below this there is 
a single extremal state at m = n = Na + Nb — 2 (not present in the 2 x Nb systems) with a complete 
set of product vectors in Ker p and no product vector in Im p. For even lower ranks the states are also here 
exceptional and do not fit into this picture. 

When we exclude these lowest-rank states, the numbers we find in the lists are in fact identical to 
the numbers of product vectors expected for generic subspaces of the given dimension. To show this we 
consider the conditions for a product state <j> <£> \ to be present in a randomly chosen subspace of dimension 
d in the Hilbert space T~L — Ha €5 T-Lb- The product state will satisfy a set of constraint equations of the 
form 

^(0®X)=Q, fc = l,2,...,JV-d (15) 

where tpk is a linearly independent set of states that are all orthogonal to the chosen d-dimensional space. 
The solutions to the equation will generally depend on a number of continuous parameters that can be 
determined by parameter counting. Thus, the product state will be specified by Na + Nb — 2 complex 
parameters (when the complex normalization factors are not included), and this number is reduced by the 
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Figure 3: Diagrams for the ranks of PPT states in the composite systems of dimensions 2x4 and 2x5. For explanations 
of the diagrams see Figs. 1 and 2. 



N — d complex constraint equations to give for the solution a remaining set of p complex parameters, with 

p = N A + N B -2- N + d (16) 
We may then distinguish between three cases, 

1) p > 0, which means d > N — Na — Nb + 2: The set of equations (15) is underdetermined, and there is 
an infinite set of product vectors in the d-dimensional subspace, described by p complex free parameters. 

2) p = 0, which means d = N — Na — Nb + 2: The number of equations matches the number of parameters 
to give a finite set of solutions. 

3) p < 0, which means d < N — Na — Nb + 2: The set of equations is overdetermined and there is in the 
generic case no solution. There may however be solutions for specially selected subspaces. 

For the case 2) the number of product vectors is given by the expression 

_ f N A + N B -2\ _ (Na+Nb- 2)1 
ps ~\ N A -1 J ~ (N A - 1)\(N B - 1)! 

The problem of finding this number of product vectors can be related to the problem in algebraic geometry 
of finding the degree of the variety defined by the Segre embedding between projective spaces P n a-"l x 
pN B -i _^ pN A N B -i [13]_ xhe number given above is identical to this degree. 

It is straight forward to check that all the numbers in the tables are consistent with these results, except 
for the special low-rank states. Apart from these states, the PPT states we find in the numerical searches 
are therefore also in this respect typical states for the given ranks (to, n). 



3.4 The lowest-rank extremal and entangled PPT states 

For the systems with Na > 2 and Nb > 2 the lowest-rank extremal PPT states with full local ranks seem 
to have a special status. These states have symmetric ranks m — n, and for these values of m (or n) we 
find no states with asymmetric ranks. They are also special in the sense that they are the only ones with no 
product vector in Im p and a finite, complete set of product vectors in Ker p. 

For the 3x3 system these are rank (4, 4) states, and the presence of product vectors in Kerp but not in 
Im p indicates that they are related to a special construction of low rank entangled PPT states for this sys- 
tem, with the use of unextendible product bases (UPB) [14]. The UPB is a set of orthogonal product vectors 
that spans a subspace (Ker p), and that cannot be extended with additional orthogonal product vectors (in 
Im p). The (4, 4) states that we find by our method do not directly exemplify this construction, since the 
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Figure 4: Diagrams for the ranks of PPT states in the composite systems of dimensions 2x2 and 2x3. For these systems 
all PPT states are separable, but the distinction between the special set of states with low, symmetric ranks m — n and 
the set of states with higher, asymmetric ranks is similar to what is found in the higher-dimensional systems. 



product vectors in Ker p are non-orthogonal, but there is a general connection to the UPB construction that 
we discuss in detail in a separate publication [11]. 

We find similar types of extremal PPT states in all the systems we have studied with dimensions of the 
subsystems larger than 2. This we take as an indication for the presence of such lowest-rank extremal states 
in all higher-dimensional systems. When we combine the assumption of a finite, complete set of product 
vectors in Ker p and no product vector in Im p with the expectations for the number of product vectors 
in generic subspaces of a given dimension, we find that these states should generally have ranks equal to 
m = n = Na + Nb — 2. We will phrase the assumption about these states in the form of the following 
conjecture: 

The lowest-rank extremal PPT state with full local ranks in an Na X Nb system, with Na and Nb 
larger than 2, is characterized by symmetric ranks m=n for the density matrix and its partial transpose, 
with value m = n = Na + Nb — 2. 

We write this as the following inequality for the rank of extremal PPT states with full local ranks, 

m, n > N A + N B - 2 (conjecture) (18) 
and compare it to the lower bound for entangled PPT states with full local ranks, 

m,n> max{iV A , N B } + 1 (HLVC) (19) 

For systems of dimension 2 x Nb we note that the HLVC bound lies above the lower bound for extremity, 
and therefore the bound (18) cannot be saturated. That is in accordance with the lack of the special type of 
lowest-rank extremal PPT state for these systems. For systems of dimension 3 x Nb (with Nb > 3) the 
two bounds coincide, and the lowest rank extremal states that we have found for these systems are indeed 
also the lowest rank entangled states. For systems of dimension 4 x Nb (with Nb > 4) there is a difference 
of 1 between the two lower bounds with the HLVC bound being the lowest. We have examined only one 
example of such systems, namely the 4x4 system. For this system we do find entangled PPT states with 
lower rank than the lowest-rank extremal PPT state. For general systems Na x Nb, with Na < Nb, the 
difference between the two bounds is according to our conjecture equal to Na — 3, and therefore increases 
linearly with the lowest dimension of the two subsystems. 

The entangled PPT states with lower ranks than the lowest rank extremal PPT states with full local ranks 
will include only extremal PPT states with less than full local ranks, when written as convex combinations 
of extremal states. We have found examples of such states with ranks (5, 5) in the 4x4 system. In the table 
such a state p appears with dim T — 2, which means that it is a convex combination of two extremal PPT 
states, 

p = (1 - x)p e + xp p (20) 

with < x < 1. One of these is a pure product state p p = wvr , where w = u ® v is the single product 
vector found in Im p. The other one is a rank (4, 4) state p e , which we identify as an extremal PPT state 
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of a 3x3 subsystem. We find this state by subtracting p, p with the value of the coefficient x determined as 
discussed in [12]. To summarize, the entangled rank (5, 5) PPT state p has full local ranks and saturates 
the HLVC bound, but is a convex combination of two extremal PPT states that both have less than full local 
ranks. 

A closer look at the above decomposition of the (5, 5) state p, motivates the following general con- 
struction of states that saturate the HLVC bound in successively higher dimensions. Assume that 

n A =Ui®u 2 , n B = v 1 ®v 2 (2i) 

where dim^i = dimVi = 3 and dimW 2 = dimV 2 = 1. We assume that U\ andU 2 are complementary 
but not necessarily orthogonal subspaces of Ha, and that Vi and V 2 are complementary but not necessarily 
orthogonal in Hb- 

Let p e be a rank (4, 4) extremal PPT state on the 3 x 3 dimensional subspace U\ <g> Vi C H. A property 
of p e is that there are no product vectors in Im p e C U\ ® Vi . Let p p be a pure product state, p p = ww^ 
with w = u®v,u ei/ 2 and v € V 2 . Then define p = (1 — x)p e + xp p with < x < 1. The image of p is 

Imp = Imp e lmp p = Im p e © (U 2 O V 2 ) C (Wi <g> Vi) (W 2 ® V 2 ) (22) 

Since U 2 ® V 2 is one dimensional, it follows that the rank of p is one higher than the rank of p e . Similarly, 
the rank of p p is one higher than the rank of pf. 

The only product vector in Im p is now w = u ® v. To see this, consider the most general product 
vector 

w' = (ui + u 2 ) ® (vi + v 2 ) = (tii <S> t>i) + (mi ® w 2 ) + (u 2 ® ui) + (u 2 (?) w 2 ) (23) 

with Ui e and e Vj. It is a sum of four vectors belonging to four complementary subspaces 
Ui ® Vj C H with i, j = 1, 2. In order to have w' e Imp we must have 

vi € (Wi®Vi)©(W 2 ®V 2 ) (24) 

but this requires that Ui ® u 2 = u 2 ® «i = 0. Since we want to have w' ^ 0, the only possibilities left are 
w' = u\ (g) «i, which is not in Im p, or w' = u 2 ® f 2 , which is not new. 

It should be clear from this analysis that the same construction can be applied to higher dimensional 
systems. For example, we may reinterpret p e to be the (5, 5) state p of the 4x4 system, as constructed 
above. We reinterpret p p to be a new pure product state that increases the local ranks from 4 to 5, in a 
similar manner as discussed, so that p = (1 — x)p e + xp p is a (6, 6) state of the 5x5 system. The analysis 
of product vectors in Im p can be done in the same way, with one exception. In this case there is already 
one product vector in Im p e and therefore there are two possible solutions for the product vector w' , namely 
the newly added vector w = u<E> v and the previous product vector in Im p e . Since the number of product 
vectors is lower than dim Im p the state is entangled, and it saturates the HLVC bound. 

The construction can be continued to arbitrarily high dimension, and it will always create an entangled 
state that saturates the HLVC bound. The state will be a convex combination of extremal PPT states 
with less than full local ranks, but it will itself have full local ranks. The construction is not restricted to 
symmetric cases, Na — Nb, since asymmetric systems can be reached by introducing product vectors 
which increase the dimension of only Ha or Hb- In this way it is possible to saturate the HLVC bound in 
all higher-dimensional systems. 

3.5 Checking for entanglement at higher ranks 

All the states with ranks above the upper limit for extremity, i.e., with 

m 2 + n 2 > N 2 + 1 (25) 

we find to have a complete set of product vectors in their image. This means that they satisfy the necessary 
condition for separability given by the range criterion, which however, does not exclude the possibility that 
they may be entangled. To get some further information about this we have made use of a separability 
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criterion introduced in [12]. This condition for separability can be seen as a strengthened version of the 
range criterion, and is based on a relation that, for separable states, exists between product vectors in Im p 
and in Imp p . We describe below the basis for this criterion and further describe the method we have 
applied for checking the criterion. 

Assume p to be a separable density operator, which therefore can be written as a convex combination 
of product states, 

P = ^Pk^ki>l (26) 
k 

with ipk — 4>k®Xk- The partially transposed density operator can then be written as 

P P = Y / P^l (27) 

k 

with tpf. = <j>t~ ® xl> where \% is the complex conjugate of \k with respect to the same basis in Hb that 
is used for the partial transposition. Therefore, corresponding to the set of product vectors {ipk}, which 
spans Im p, there is a set of product vectors {ipk}, which we will refer to as the conjugate set, that spans 
lmp p . This implies that a necessary condition for separability of a density operator p is that the number 
of pairs of product vectors in Im p and in Im p p that are conjugate, is equal to or larger than the ranks of 
both p and p p . We write the condition as 

K > max{m, n} (28) 

with K as the number of conjugate pairs of product vectors in Im p and Im p p . 

The above condition for separability is effectively restricted to cases where the ranks m and n of p and 
p p are not too large. This follows since, if the dimensions of Imp and Imp p are sufficiently large, the 
number of such pairs will necessarily be infinite. The condition for the number of pairs to be finite can be 
determined by essentially the same method of counting parameters as used to determine the typical number 
of product vectors in a Hilbert space of given dimension, as discussed in Sect. 3.3. Thus, a conjugate pair 
of product vectors, (ip, tp), with ip = <p ® \ an d V' — 4> ® X*> has to satisfy two sets of equations, 

0^ = 0, g$ = (29) 

where {8i} is a basis of Ker p and {^} is a basis of Ker p p . The number of equations that ip has to satisfy 
will be equal to, or will exceed the number of free parameters in this product state, provided the following 
condition is satisfied [12] 

dim Ker p + dim Ker p p > dim "Ha + dim"Hs - 2 (30) 

If the condition is satisfied with proper inequality there will typically be no solution and if it is satisfied 
with equality there will be a finite set of solutions. Expressed in terms of the ranks m and n and the Hilbert 
space dimensions 7Y, Na and Nb the inequality takes the form 

m + n< 2N-N A -N B + 2 (31) 

and this sets the limit for applications of the criterion. In the diagrams this limit is indicated by the dashed 
green line. 

For states that satisfy the above inequality, we have found a practical method to check the separability 
condition (28) by applying essentially the same double-eigenvalue method as used for detecting product 
vectors in Im p (and in Ker p), and which is described in Appendix A. The outline of the method is the 
following. Let P be the orthogonal projection on Imp and Q the orthogonal projection of Imp p . We 
consider the following bilinear function of the product state ip = <j> ® x> 

/(V) = ft (i - P)i> + ft (i - Q)i> (32) 

and search for the minima of the function where / = 0. Since both terms in the expression (32) are non- 
negative, such a solution will give zero for each term separately, and from this follow the relations Pip = tp 
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and Qip = ip, which are equivalent to the two sets of equations (29). The second term in (32) can be 
rewritten in the following way, — Q)ip = tp^(l — Q p )tp, with Q p as the partial transpose of Q. This 
gives 

/(V>) =V f (21- P-Q P )^ (33) 

and written in this way the function f{ifj) has the same form as the function that is minimized in the search 
for product vectors in Imp (see Sect. 2.3). The only difference is that the operator 1 — P is replaced by 
21 — P — Q p . The same method to search for product vectors with vanishing value for the function can 
therefore be used. 

The result is that for all the states we have found with m > Na + Nb — 1, and where (31) is satisfied as 
a proper inequality, there are no pairs of conjugate product vectors, and the states are therefore entangled. 
In the diagrams these states are located below the dashed green lines. For the states on the dashed green 
line in the diagrams the condition (31) is satisfied with equality, and for these states we find a finite number 
K of pairs of conjugate product vectors that is sufficiently large to satisfy the separability condition (28). 
This means that just counting the number of pairs of conjugate product vectors is insufficient to determine 
if the states are entangled. However, with only a finite number of product vectors available it is possible to 
check whether the density operator can be reconstructed as a convex combination of these product states. 
The result is that for all ranks where (31) is satisfied with equality we find entangled states. However in a 
small number of cases we find both separable and entangled states with the same set of ranks (to, n). That 
happens for the (6, 6), (5, 7) and (7, 5) states of the 2x4 system. 

4 Concluding remarks 

The results presented in this paper are based on the use of a numerical method to search for density matrices 
p of a composite, bipartite system with specified values for the ranks of p and its partial transpose p p . The 
method works well for systems where the Hilbert space dimension N is not too large, in our calculations 
with N < 20. In higher dimensions the main problem is a too slow convergence of the iteration procedure. 
Additional methods to find the number of product vectors in the image and kernel of the density matrices 
and to determine the dimension of the corresponding face of the set of PPT states have been described. The 
results based on the use of these methods are listed in the tables and displayed in the diagrams. 

The results obtained for the low-dimensional systems that we have studied reveal several regularities. 
For sufficiently low ranks we find only separable states of a specific form. Above a certain value for the 
ranks of p and p p the states we find are typically extremal PPT states, until the ranks reach an upper limit. 
In our discussion we suggest that there are in fact two lower bounds, one for entangled states with full local 
ranks and another, generally more restrictive, for extremal PPT states with full local ranks. The first one we 
identify as the bound on entangled PPT states discussed in [12], and we show, by an explicit construction, 
how this bound can be saturated. Based on certain properties of the extremal states with minimal ranks that 
are found in our searches, we conjecture a specific value for the second one, the lower bound on the ranks 
of extremal PPT states. The property we focus on is the number of product vectors in the image and in 
the kernel of this state, which makes these states different from the higher rank extremal states. Assuming 
this property to be present for the lowest rank extremal states in general, we draw the conclusion about the 
lower bound. 

Above this lower bound we find in our searches a large set of ranks (to, n) of p and p p where the states 
are extremal. There is an upper limit to the ranks of these states, which can be understood as following from 
a geometrical constraint on the face of the set V to which an extremal PPT state belongs. It is of interest to 
note that all the states we find for sufficiently low ranks are limited to the symmetric cases m — n, whereas 
for to and n above this limit we find in addition states for essentially all the asymmetric values of the ranks. 
However, as we have stressed, the states we find by our method should be considered as typical states for 
the given ranks (to, n). This means that we cannot exclude the presence of untypical states also for ranks 
where we have not identified any PPT state. 

Concerning this last point it is of interest to relate the results for the properties of the PPT states 
found in our searches with those of other entangled and extremal PPT states referred to in the literature 
[19, 14, 15, 16, 17, 18]. These states are based on special constructions which lead to certain classes of 
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extremal PPT states. Among the states presented by these special constructions we have found no example 
of states with ranks different from those referred to in our tables and diagrams. The main difference, 
however, is that these specially constructed states are not necessarily typical in the meaning used here, and 
in particular the number of product vectors in the image and kernel may be larger than the minimal values 
that we find in our searches. 

In the discussion of our results we have put some emphasis on the properties of the lowest rank extremal 
PPT states. For all the systems we have studied they are special in the sense that they have no product state 
in their image, but a complete, finite set of product vectors in their kernel. In a separate publication [11] 
we have made a detailed study of these states for the 3x3 system, where they have ranks (4, 4). We show 
there that these states can be related to states constructed from unextendible product bases, and the set of 
such states can be given an explicit parameterization. We do not know how to generalize this construction 
to higher dimensions, and the question of a general parameterization of extremal PPT states remains as an 
interesting problem for future work. 
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A The minimum double eigenvalue problem 

Given the Hermitian matrix A, the problem considered here is to minimize the expectation value tp^Atp 
over product vectors = <fi ® x> with tne normalization conditions <ft<f> = x*X = 1- We present an 
iteration method which may converge to different local minima, depending on the starting point for the 
iterations. We introduce Lagrange multipliers A, [i and define 



i,j,k,l \ i / \ 3 ) 

The following equations must hold at the minimum, or more generally at an extremum, 



(34) 



= XjAij;ki<f>kXi ~ A0 4 = , — = <f>* A ij;ki<PkXi ~ MX; = , (35) 
j,k,l X i i,k,l 



with 



Y <t>*<f>i = Y X*3*i = 1 ' A = V = Y fiiXjAijw&kXi • (36) 

i j i,j,k,l 

Given an approximate solution ip = <j> ® \> we compute a next approximation ip' = <f>' ® \' as follows. 
We compute x = Bcj) — A</>, y = C\ — Ax, with A = ip^Aip and 

Bik = Y X*jAij-,kiXi i Cji = Y ^*i A ^-M<i>k ■ (37) 

j,l i.k 

In practice, we compute z = Atp, Ui = J2jXj z ij> v j = ^i^i^i' ^ = ^ u — X* v i x ~ u — \<p, 
y = v - Xx- Then <\>' = Ni (4> + ex), \' = N 2 (x + ty), where Ni , N 2 are normalization factors, and where 
e is determined in the following way. To first order in e we have 

((f> + ex) ®(x + ey) = ip + ew , (38) 

with w = 4>®y + x®x- Note that <p^x = x^V — 0, and hence ip^w = 0. We compute e by minimizing 
s^As/s^s with s = ip + ew. This is an eigenvalue problem in the two dimensional subspace spanned by ip 
and w, and it can be solved analytically. 



13 



B Tables 



2x4 



(m,n) 


m 2 + n 2 - 


N 2 


dimJF 


(rA,r B ) 


#pv [Im p] 


#pv [Ker p] 


(8,8) 


64 




64 


(2,4) 


oo/8 





(8 7) 


49 




49 


(2 4) 


oo/8 





(8,6) 


36 




36 


(2,4) 


oo/8 





(7,7) 


34 




34 


(2,4) 


oo/7 





(8,5) 


25 




25 


(2,4) 


oo/8 





(7,6) 


21 




21 


(2,4) 


oo/7 





(7,5) 


10 




10 


(2,4) 


oo/7 





(6,6) 


8 




8 


(2,4) 


oo/6 





(6,5) 


-3 




1 


(2,4) 


oo/6 





(5,5) 


-14 




1 


(2,4) 


oo/5 





(4,4) 


-32 




4 


(2,4) 


4/4 


4/4 


(3,3) 


-46 




3 


(2,3) 


3/3 


oo/5 


(2,2) 


-56 




2 


(2,2) 


2/2 


oo/6 


(1,1) 


-62 




1 


(1,1) 


1/1 


oo/8 


2x5 


(m,n) 


m 2 + n 2 - 


N 2 


dimJF 


(r A ,r B ) 


#pv [Im p] 


#pv [Ker p] 


(10,10) 


100 




100 


(2,5) 


oo/lO 





(10,9) 


81 




81 


(2,5) 


oo/lO 





(10,8) 


64 




64 


(2,5) 


oo/lO 





(9,9) 


62 




62 


(2,5) 


oo/9 





(10,7) 


49 




49 


(2,5) 


oo/lO 





(9,8) 


45 




45 


(2,5) 


oo/9 





(9,7) 


30 




30 


(2,5) 


oo/9 





(8,8) 


28 




28 


(2,5) 


oo/8 





(9,6) 


17 




17 


(2,5) 


oo/9 





(8,7) 


13 




13 


(2,5) 


oo/8 





(8,6) 







1 


(2,5) 


oo/8 





(7,7) 


-2 




1 


(2,5) 


oo/7 





(7,6) 


-15 




1 


(2,5) 


oo/7 





(6,6) 


-28 




1 


(2,5) 


oo/6 





(5,5) 


-50 




5 


(2,5) 


5/5 


5/5 


(4,4) 


-68 




4 


(2,4) 


4/4 


oo/6 


(3,3) 


-82 




3 


(2,3) 


3/3 


oo/7 


(2,2) 


-92 




2 


(2,2) 


2/2 


oo/8 


(1,1) 


-98 




1 


(1,1) 


1/1 


oo/9 



Table 1: Numerical results for the 2x4 and the 2x5 systems. The first column lists the ranks of p and p p 
where PPT states have been found. The second column lists the lower limit for the value of the dimension of 
the face of JC(V) for the given ranks (to, n), while the third column lists the actual values of the dimensions 
for the states we have found. The fourth column lists the values of the local ranks with respect to subsystems 
A and B. The fifth and sixth columns give the number of product vectors in Imp and Kerp, respectively. 
In each of these columns two numbers are given, with the number to the left as the total number and the 
one to the right as the number of linearly independent product vectors, oo indicates that we find no upper 
limit to the number of product vectors that can be generated. In the present tables the extremal PPT states 
with no product state in Im p and a complete set in Ker p, which we find in all the other tables, are missing. 
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3x3 



( m n ^ 


m 2 + n 2 — 


N 2 


dim 


( v a r-D 'l 


#pv [Im p] 


#pv [Ker p] 


(9,9) 


81 




81 


(3,3) 


oo/9 





(9,8) 


64 




64 


(3,3) 


oo/9 


u 


(9,7) 


/in 

49 




a a 
49 


(3,3) 


oo/9 


U 


(8,8) 


An 
4 / 




An 
4 / 


(3,3) 


oo/8 


u 


(9,6) 


36 




36 


(3,3) 


oo/9 





(o,7) 


32 




32 


(3,3) 


oo/8 


U 


(8,6) 


1 a 

iy 




1 a 
19 


(3,3) 


oo/8 


u 


(1,1) 


1 n 
1 / 




1 n 
1 / 


(3,3) 


In 

OO// 


u 


(8,5) 


8 




8 


(3,3) 


oo/8 





(7,6) 


4 




4 


(3,3) 


oo/7 





(7,5) 


-7 




t 
1 


(3,3) 


In 

OO// 


A 
U 


(6,6) 


-9 




1 


(3,3) 


oo/6 





(6,5) 


-20 




1 


(3,3) 


oo/6 





(5,5) 


-31 




1 


(3,3) 


6/5 





(4,4) 


-49 




1 


(3,3) 





6/5 


(3,3) 


-63 




3 


(3,3) 


3/3 


oo/6 


(2,2) 


-73 




2 


(2,2) 


2/2 


oo/7 


(1,1) 


-79 




1 


(1,1) 


1/1 


oo/8 


3x4 


(m,n) 


> > 
m 2 + n 2 - 


tvtV 

N 2 


dim T 


(r A ,r B ) 


#pv [Im p] 


#pv [Ker p] 


12,11) 


121 




121 


(3,4) 


oo/12 





(12,10) 


100 




100 


(3,4) 


oo/12 





(11,11) 


98 




98 


(3,4) 


oo/ll 





(12,9) 


81 




81 


(3,4) 


oo/12 





(11,10) 


77 




77 


(3,4) 


oo/ll 





(12,8) 


64 




64 


(3,4) 


oo/12 





(11,9) 


58 




58 


(3,4) 


oo/ll 





(10,10) 


56 




56 


(3,4) 


oo/lO 





(11,8) 


41 




41 


(3,4) 


oo/ll 





(10,9) 


37 




37 


(3,4) 


oo/lO 





(11,7) 


26 




26 


(3,4) 


oo/ll 





(10,8) 


20 




20 


(3,4) 


oo/lO 





(9,9) 


18 




18 


(3,4) 


oo/9 





(11,6) 


13 




13 


(3,4) 


oo/ll 





(10,7) 


5 






(3,4) 


oo/lO 





(9,8) 


1 






(3,4) 


oo/9 





(10,6) 


-8 






(3,4) 


oo/lO 





(9,7) 


-14 






(3,4) 


oo/9 





(8,8) 


-16 






(3,4) 


oo/8 





(9,6) 


-27 






(3,4) 


oo/9 





(8,7) 


-31 






(3,4) 


oo/8 





(8,6) 


-44 






(3,4) 


oo/8 





(7,7) 


-46 






(3,4) 


10/7 





(7,6) 


-59 






(3,4) 


10/7 





(6,6) 


-72 






(3,4) 








(5,5) 


-94 






(3,4) 





10/7 


(4,4) 


-112 




4 


(3,3) 
(3,4) 




4/4 


oo/8 
oo/8 


(3,3) 


-126 




3 


(3,3) 


3/3 


oo/9 


(2,2) 


- 136 




2 


(2,2) 


2/2 


oo/lO 


(1,1) 


-142 




1 


(1,1) 


1/1 


oo/ll 



Table 2: Numerical results for the 3x3 and the 3x4 systems. For explanations see Table 1. Except for one 
case in the 3x4 system we find only one type of state, characterized by the listed properties, for each set of 
ranks (m, n). The exception is the (4, 4) where we find both extremal PPT states with less than full local 
ranks and separable states of the same construction a those of lower ranks. 
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3x5 



(m, n) 


m 2 + n 2 - N 2 


dim T 


(fa, i"b) 


#pv [Im p] 


#pv [Ker p] 


f 1 5 14) 


196 


196 


n 5) 


oo/15 





(15,13) 


169 


169 


n 5) 


oo/15 





H4 141 


167 


167 


n 5) 


oo/14 





(\ 5 12) 


144 


144 


n 5) 


oo/15 





(14 13") 


140 


140 


n 5) 


oo/14 





(15 11) 


121 


121 


n 5) 


oo/15 





CI 4 12) 


115 


115 


n 5) 


oo/14 





(13,13) 


113 


113 


n 5) 


oo/13 





C14 11) 


92 


92 


n 5) 


oo/14 





en 12) 


88 


88 


n 5) 


oo/13 





(14 10) 


71 


71 


n 5) 


oo/14 





(13 11) 


65 


65 


n 5) 


oo/13 





(12,12) 


63 


63 


n 5) 


oo/12 





(14 9) 


52 


52 


n 5) 


oo/14 





(13,10) 


44 


44 


n 5) 


oo/13 





C12 11) 

v iz -, 1 A / 


40 


40 


n 5) 


oo/12 





("14 8) 


35 


35 


n 5) 


oo/14 





en 9) 


25 


25 


n 5) 


oo/13 





CI4 7) 


20 


20 


n 5) 


oo/14 





(12,10) 


19 


19 


n 5) 


oo/12 





n 1 id 


17 


17 


n 5) 


oo/ll 





en 8) 


8 




n 5) 


oo/13 





("12 9) 





j 


n 5) 


oo/12 





C1 1 10) 


_4 


j 


n 5) 


oo/ll 





(13 7) 


_7 


j 


n 5) 


oo/13 





(12 8) 


-17 


j 


n 5) 


oo/12 





C1 1 9) 

V A A , 7 / 


-23 


j 


n 5) 


oo/ll 





(10,10) 


-25 


j 


n 5) 


oo/lO 





(12 7) 


-32 


j 


n 5) 


oo/12 





(1 1 8) 

V A A »°' 


-40 


j 


n 5) 


oo/ll 





no 9) 


-44 


j 


n 5) 


oo/lO 





(11,7) 


-55 


1 


(3,5) 


oo/ll 





no 8) 


-61 




(3 5) 


co/10 





(9 9) 


-63 


-. 


C3 5) 


15/9 





no 7) 


-76 


j 


(3 5) 


co/10 





(9 8) 


-80 


j 


(3 5) 


15/9 





(9,7) 


-95 


j 


(3,5) 


15/9 





(8,8) 


-97 




(3,5) 








(8,7) 


-112 




(3,5) 








(7,7) 


-127 




(3,5) 








(6,6) 


-153 




(3,5) 





15/9 


(5,5) 


-175 


5 


(3,4) 
(3,5) 




5/5 


00/10 
00/10 


(4,4) 


-193 


4 


(3,4) 


4/4 


00/11 


(3,3) 


-207 


3 


(3,3) 


3/3 


00/12 


(2,2) 


-217 


2 


(2,2) 


2/2 


00/13 


(1,1) 


-223 


1 


(1,1) 


1/1 


00/14 



Table 3: Numerical results for the 3x5 systems. For explanations see Table 1. We here find two types of 
states with ranks (5, 5). 
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4x4 



(m, n) 


in + n — IN 


dim 7~ 


(TA,rB) 


#pv [Im p] 


#pv [Ker p] 


(16,16) 


256 


256 


(4,4) 


oc/16 





(16,15) 


225 


225 


(4,4) 


oc/16 





(16,14) 


196 


196 


(4,4) 


oc/16 





(15,15) 


194 


194 


(4,4) 


oc/15 





(16,13) 


169 


169 


(4,4) 


oc/16 





(15,14) 


165 


165 


(4,4) 


oc/15 





(16,12) 


144 


144 


(4,4) 


oc/16 





(15,13) 


138 


138 


(4,4) 


oc/15 





(14,14) 


136 


136 


(4,4) 


oc/14 





(16,11) 


121 


121 


(4,4) 


oc/16 





(15,12) 


113 


113 


(4,4) 


oc/15 





(14,13) 


109 


109 


(4,4) 


oc/14 





(15,11) 


90 


90 


(4,4) 


oo/15 





(14,12) 


84 


84 


(4,4) 


oc/14 





(13,13) 


82 


82 


(4,4) 


oc/13 





(15,10) 


69 


69 


(4,4) 


oc/15 





(14,11) 


61 


61 


(4,4) 


oc/14 





(13,12) 


57 


57 


(4,4) 


oc/13 





(15,9) 


50 


50 


(4,4) 


oo/15 





(14,10) 


40 


40 


(4,4) 


oc/14 





(15,8) 


33 


33 


(4,4) 


oc/15 





(13,11) 


34 


34 


(4,4) 


oc/13 





(12,12) 


32 


32 


(4,4) 


oc/12 





(14,9) 


21 


21 


(4,4) 


oc/14 





(15,7) 


18 


18 


(4,4) 


oc/15 





(13,10) 


13 


13 


(4,4) 


oc/13 





(12,11) 


9 




(4,4) 


oc/12 





(14,8) 


4 




(4,4) 


oc/14 





(13,9) 


-6 


1 


(4,4) 


oc/13 





(14J) 


-11 


1 


(4,4) 


oc/14 





(12,10) 


-12 


1 


(4,4) 


oo/12 





(11,11) 


-14 


1 


(4,4) 


oc/11 





(13,8) 


-23 


1 


(4,4) 


oc/13 





(12,9) 


-31 


1 


(4,4) 


oc/12 





(11,10) 


-35 


1 


(4,4) 


oo/ll 





(13,7) 


-38 


1 


(4,4) 


oc/13 





(12,8) 


-48 


1 


(4,4) 


oc/12 





(11,9) 


-54 


1 


(4,4) 


oo/ll 





(10,10) 


-56 


1 


(4,4) 


20/10 





(12,7) 


-63 


1 


(4,4) 


oc/12 





(11,8) 


-71 


1 


(4,4) 


oo/ll 





(10,9) 


-75 


1 


(4,4) 


20/10 





(11,7) 


-86 


1 


(4,4) 


oo/ll 





(10,8) 


-92 


1 


(4,4) 


20/10 





(9,9) 


-94 




(4,4) 








(10,7) 


-107 


1 


(4,4) 


20/10 





(9,8) 


-111 




(4,4) 








(9,7) 


-126 




(4,4) 








(8,8) 


-128 


1 


(4,4) 








(8,7) 


-143 


1 


(4,4) 








(7,7) 


-158 


1 


(4,4) 








(6,6) 


-184 


1 


(4,4) 





20/10 


(5,5) 


-206 


2 
5 


(4,3) 
(4,4) 
(4,4) 




1/1 
5/5 


oo/ll 
oo/ll 
oo/ll 


(4,4) 


-224 


1 
4 


(3,3) 
(4,4) 




4/4 


0O/12 
00/12 


(3,3) 


-238 


3 


(3,3) 


3/3 


00/13 


(2,2) 


-248 


2 


(2,2) 


2/2 


00/14 


(1,1) 


-254 


1 


(1,1) 


1/1 


00/15 



Table 4: Numerical results for the 4x4 systems. For explanations see Table 1. The pattern of the listed 
properties is much like that of the other tables, but here with more low-rank states below the lowest-rank 
(to = n — 6) extremal PPT state with full local ranks. 
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